IFIC/12-31 



Global status of neutrino oscillation parameters after Neutrino-2012 

D. V. ForeroQ M. TortolaQ and J. W. F. ValkS 

AHEP Group, Instituto de Fisica Corpuscular - C.S.I.C./Universitat de Valencia 
Edificio de Institutos de Paterna, Apartado 22085, E-46071 Valencia, Spain 

Here we update the global fit of neutrino oscillations in Refs. [f|, Q| including the recent measure- 
ments of reactor antineutrino disappearance reported by the Double Chooz, Daya Bay and RENO 
experiments, together with latest MINOS and T2K appearance and disappearance results, as pre- 
sented at the Neutrino-2012 conference. We find that the preferred global fit value of #13 is quite 
large: sin 2 #13 ~ 0.025 for normal and inverted neutrino mass ordering, with #13 — now excluded 
I at more than 10a. The impact of the new 813 measurements over the other neutrino oscillation pa- 

■ rameters is discussed as well as the role of the new long-baseline neutrino data and the atmospheric 

^ , neutrino analysis in the determination of a non- maximal atmospheric angle 823- 

O- 
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D . I. INTRODUCTION 

& 

Recent measurements of 6*13 have been reported by the reactor experiments Double Chooz Daya Bay [4| and 
RENO Q . These experiments look for the disappearance of reactor antineutrinos over baselines of the order of 1 km, 
qq , due to neutrino oscillations mainly driven by the third mixing angle #13 of the lepton mixing matrix 0, Q] • Up to now 
the most sensitive measurements of reactor antineutrinos were reported by the past reactor experiments CHOOZ Q 
and Palo Verde jgj]. 

Compared to their predecessors, the new reactor experiments have larger statistics, thanks to their increased reactor 



o 



, power and the bigger antineutrino detector size. On the other hand, one of their most important features is that they 
have detectors located at different distances from the reactor core. As a result measurements at the closest detectors 
can be used in order to predict the expected event number at the more distant detectors, avoiding the need to rely 
on theoretical calculations of the produced antineutrino flux at the reactors. As a consequence these experiments 
have for the first time observed the disappearance of reactor antineutrinos over short distances, providing the first 
^ ' measurement of the mixing angle #13, so far unknown. 

Last year there were some indications for a non-zero #13 mixing angle coming from the observation of electron 
neutrino appearance on a muon neutrino beam at the accelerator oscillation experiments T2K [nl | and MINOS 11 1 . 



Together with the hints from the solar and atmospheric neutrino data samples, the global analysis of neutrino oscilla- 
tion data reported indications of non-zero #13 between 3 and 4o\ depending on the treatment of short-baseline reactor 
data in the full analysis (see Refs P, Q for more details 1 ). Here we update the global fit of neutrino oscillations given 
111 Refs. U0 by including the recent measurements of reactor antineutrino disappearance reported by the Double 
Chooz, Daya Bay and RENO experiments. In addition to the above-mentioned data our analysis includes also the 
most recent reactor neutrino data reported at the Neutrino-2012 Conference , [l^|, as well as the latest results of 
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the MINOS 16| and T2K 17,[18| experiments. The role of the new $13 measurements upon the determination of 
the remaining oscillation parameters is also analyzed. Particularly important is the impact of the new accelerator 
neutrino data as well as the details of the atmospheric neutrino analysis upon the determination of the atmospheric 
mixing angle #23- 



II. REACTOR EXPERIMENTS: DOUBLE CHOOZ, DAYA BAY AND RENO 



This year, in chronological order, the Double Chooz, Daya Bay and RENO Collaborations have reported measure- 
ments of the electron antineutrino disappearance with important levels of statistical significance. 

The Double Chooz (DC) experiment, located in France, is a reactor experiment planned to have two detectors 
and two reactors. In its first stage DC has reported 101 days of running 3|, with only the far detector operating 
so far. The near detector (ND) is expected to start operation by early 2013. The two reactors are approximately 
equal, with an individual power of 4.25 GW t h and are placed at a distance of 1050 m from the far detector. The 
detector has a fiducial volume of 10 m 3 of neutrino target liquid. From the analysis of the rate and the energy 
spectrum of the prompt positrons produced by the reactor antineutrinos, the DC collaboration find sin 2 29\ 3 = 
0.086 ± 0.041(stat) ± 0.030(syst). Using only the ratio of observed to expected events they get a slightly higher best 
fit value: sin 2 2#i3 = 0.104 ± 0.030(stat) ± 0.076(syst). A more recent analysis of DC data with an exposure of 227.93 
live days 14j has reported the observation of 8249 candidate electron antineutrinos while 8937 were expected in the 
absence of oscillations. Using a rate plus spectral shape analysis the following best fit value for the reactor angle is 
obtained: sin 2 29 13 = 0.109 ± 0.030(stat) ± 0.025(syst). 

The Daya Bay (DYB) reactor experiment 0| is a neutrino oscillation experiment designed to measure the mixing 
angle #13 as well. The experiment is placed in China and it contains an array of three groups of detectors and three 
groups of two-reactor cores. The far group of detectors (far hall) is composed of three detectors and the two near 
halls are composed by one and two detectors, respectively. In order to reduce systematic errors, the detectors are 
approximately equal, with a volume of 20 tons of Gadolinium-doped liquid scintillator as neutrino target material. 
The reactor cores are approximately equal as well, with a maximum power of 2.9 GWth (total power of 17.4 GWth) 
and the distances to the detectors range from 350 to 2000 m approximately. The rate-only analysis performed by 
the DYB collaboration finds a best fit value of sin 2 26»i 3 = 0.092 ± 0.016(stat) ± 0.005(syst). A zero value for 13 
is excluded with a significance of 5.2c. New results presented in the Neutrino-2012 Conference [3 with 2.5 times 
more statistics allow a stronger rejection for 6*13 = that now is excluded at almost 8a by DYB alone. A rate-only 
statistical analysis of the new DYB data reports a best fit value of sin 2 26*13 = 0.089 ± 0.010(stat) ± 0.005(syst). 

The RENO experiment 0| is situated in South Korea and it has been running for 229 days. It shares some features 
with DC and DYB. RENO has six reactor cores, distributed along a 1.3 km straight line. Two of the reactors have 
a maximum power of 2.66 GWth while the other four may reach 2.8 GWth- Reactor antineutrinos are detected by 
two identical detectors, labeled as near and far, located at 294 and 1383 m from the reactor array center. Each 
RENO detector contains 16 tons of Gadolinium-doped Liquid Scintillator. Based on a rate-only analysis, the RENO 
Collaboration finds sin 2 2#i3 = 0.113 ± 0.013(stat.) ± 0.019(syst.), together with a 4.9<r exclusion for #13 = 0. 



Reactor event calculation 

Reactor antineutrinos are produced by the fission of the isotopes 235 U, 239 Pu, 241 Pu and 238 U. Each fissile isotope 
contributes to the total reactor neutrino flux and fuel content with a certain fission fraction f l that can be calculated 
through a detailed simulation of the core evolution. After their production, the reactor antineutrinos are detected at 
the experiments via the inverse beta decay process, looking for a delayed coincidence between the positron annihilation 
and the neutron capture in the target material. The window of positron energy covered by the three experiments 
described above ranges from 0.7 to 12 MeV approximately. 
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For a given experiment, the total number of events expected at the ith detector coming from the rth reactor can 
be calculated as: 



N P r 
iv p r th 
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(1) 



where N p is the number of protons in the target volume, P[ h is the total reactor power, denotes the efficiency of 
the detector and (Efi S ) = J^i f ^Us * s ^ ne avera ge energy released per fission, calculated from the individual fission 
fractions /' and the energy release per fission for a given isotope / taken from Ref. 19). For the antineutrino flux 
prediction $> r (E v ) we use the recent parameterization given in Ref. (20J as well as the new normalization for reactor 
antineutrino fluxes updated in Ref. [2l|. The inverse beta decay cross section <tibd(-Ev) is taken from Ref. [22j. Finally, 
for the neutrino propagation factor P(E Ul Li r ) we use the full three- neutrino disappearance probability. The distance 
between reactor and detector Li r is also used to correct the total antineutrino flux at the detector site. In order to 
minimize the dependence upon the predicted normalization of the antineutrino spectrum, we analyze the total rate of 
expected events at the far detector/s in the presence of oscillations over the no-oscillation prediction. This way, our 
statistical analysis is free of correlations among the different reactor data samples, since the relative measurements 
do not rely on flux predictions. 



III. GLOBAL ANALYSIS 

In our global analysis of neutrino oscillation parameters we combine the recent reactor data from Double Chooz, 
Daya Bay and RENO with all the remaining relevant experiments, as follows. 

A. Solar neutrino and KamLAND data 



We include the most recent solar neutrino data from the radiochemical experiments Homestake [2c 
Gallex/GNO 24| and SAGE 25 L as well as the latest published data from Borexino [26(, the three phases of the 



Super-Kamiokande experiment 27-2^] and the three phases from the Sudbury Neutrino Experiment SNO 30, 3lj . 
For our simulation of the production and propagation of neutrinos in the Sun we consider the most recent update of 

32j, fixing our calculations to the low metallicity model labeled as AGSS09. The impact of 



the standard solar model 

the choice of a particular solar model over the neutrino oscillation analysis has been discussed in the arXiv updated 
version of Ref. 0. We also include the most recent results published by the KamLAND reactor experiment with a 
total livetime of 2135 days, including the data collected during the radiopurity upgrade in the detector 



B. Atmospheric and accelerator neutrino data 



In our global fit we use the atmospheric neutrino analysis done by the Super-Kamiokande Collaboration [35j. The 
oscillation analysis has been performed within the one-mass scale approximation, neglecting the effect of the solar mass 
splitting and includes the atmospheric results from the three phases of the Super-Kamiokande experiment. Concerning 
the long-baseline data, we include the most recent results from the MINOS and T2K long-baseline experiments 
released last June at the Neutrino 2012 Conference. We consider the appearance and disappearance channels for both 
experiments as well as the neutrino and antineutrino data for MINOS. These new long-baseline results imply some 
improvements with respect to the previous MINOS and T2K data in Refs. [l3,[ll|,EBS|- 

On the one hand, the new 

results on — > v e appearance searches allow a better determination of the #13 mixing angle, although its current 
determination is fully dominated by the Daya Bay reactor data. On the other hand, and here lies the most relevant 
implication of the new long-baseline data, they show a preference for a non-maximal atmospheric mixing angle #23 hi 
the and channels. The impact of this preference on the determination of #23 in our global fit will be discussed 
in the next section. 
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FIG. 1: Upper panels: A^ 2 as a function of sin 2 613 from the analysis of the total event rate in Daya Bay (solid magenta/light 
line), RENO (dotted line) and Double Chooz (dashed magenta/light line) as well as from the analysis of long-baseline (dashed 
blue/dark line) and global neutrino data (solid blue/dark line). Except for the case of the global fit here we have fixed the 
remaining oscillation parameters to their best fit values. Lower panels: contours of A\ 2 = 1, 4, 9 in the sin 2 #13 — 5 plane from 
the global fit to the data. We minimize over all undisplayed oscillation parameters. Left (right) panels are for normal (inverted) 
neutrino mass hierarchy. 



C. Global fit results 

Here we summarize the results for the neutrino oscillation parameters obtained in our present global analysis. For 
details on the numerical analysis of all the neutrino samples see Refs. 0,0,0,0 and references therein. 

The results obtained for sin 2 #13 and 5 are summarized in Fig. Q] In the upper panels we show the A\ 2 profile as a 
function of sin 2 #13 for normal (left panel) and inverted (right panel) neutrino mass hierarchies. The solid blue/dark 
line corresponds to the result obtained from the combination of all the data samples while the others correspond 
to the individual reactor data samples and the combination of the long-baseline MINOS and T2K appearance and 
disappearance data, as indicated. One sees from the constraints on sin 2 #13 coming from the different data samples 
separately that, as expected, the global constraint on 13 is dominated by the recent Daya Bay measurements. For 
both neutrino mass hierarchies we find that the 3<r indication for #13 > obtained in our previous work Q due 
mainly to the first indications observed by MINOS and T2K is now overwhelmingly confirmed as a result of the recent 
reactor data. Thus, in our global fit we obtain a A% 2 ~ 104, resulting in a 10.2er exclusion of #13 = for both mass 
hierarchies. In the lower panels of Fig. [T]we show the contours of A^ 2 = 1,4,9 in the sin 2 6*13 — 5 plane from the 




FIG. 2: A\ 2 profiles as a function of all the neutrino oscillation parameters sin 2 #12, sin 2 #23, sin 2 #13, Am^, Am^ and 3. For 
the central and right panels the solid lines correspond to the case of normal mass hierarchy while the dashed lines correspond 
to the results for the inverted mass hierarchy. 



global fit to the neutrino oscillation data. In this plane we find the following best fit points: 

sin 2 6>i3 = 0.0246 , 8 = 0.80tt (normal hierarchy), (2) 
sin 2 6»i3 = 0.0250 , 8 = -0.03tt (inverted hierarchy). (3) 

In our previous analysis Q there was a "preferred region" at A% 2 = 1 for the CP phase 8 for normal neutrino mass 
ordering, as a result of the complementarity between MINOS and T2K appearance data. One sees that this effect has 
been diluted after the combination with the new reactor data, so no "preferred region" for the CP phase 8 remains at 
A% 2 = 1 2 . For this reason we marginalize over the CP phase 8 (and all other oscillation parameters), obtaining for 
the best fit, one-sigma errors, and the significance for 6*13 > 0: 

sin 2 9i3 = 0.0246+°;™ 2 ^ , A X 2 = 103.5 (10.2c-) (normal), 

sin 2 6>i 3 = 0.0250tg;E$2? , A X 2 = 104.7 (10. 2a) (inverted). ^' 

When compared with the previous analysis in Refs. we remark that here we are not including previous short 

baseline reactor experiments, which would lead to a somewhat less significant result for the exclusion of #13 = 0. 



2 Note that, given the approximations adopted in the atmospheric neutrino analysis in Ref. |3"j |. the sensitivity to the parameter S in our 
global fit comes only from long-baseline neutrino data. 
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parameter 


best fit 


la range 


2a range 


3a range 


Am|i [10" 5 eV 2 ] 


7.62 


7.43-7.81 


7.27-8.01 


7.12-8.20 


|Am§i| [10" 3 eV 2 ] 


2.55 
2.43 


2.46 - 2.61 
2.37 - 2.50 


2.38 - 2.68 
2.29 - 2.58 


2.31 - 2.74 
2.21 - 2.64 


sin 2 812 


0.320 


0.303-0.336 


0.29-0.35 


0.27-0.37 


sin 2 6 23 


0.613 (0.427) a 
0.600 


0.400-0.461 & 0.573-0.635 
0.569-0.626 


0.38-0.66 
0.39-0.65 


0.36-0.68 
0.37-0.67 


sin 2 0i3 


0.0246 
0.0250 


0.0218-0.0275 
0.0223-0.0276 


0.019-0.030 
0.020-0.030 


0.017-0.033 


S 


0.807T 
-0.037T 


- 2tt 


- 2tt 


- 2tt 



"This is a local minimum in the first octant of 823 with A\ 2 = 0.02 with respect to the global minimum 

TABLE I: Neutrino oscillation parameters summary. For Am 2 !, sin 2 823, sin 2 #13, and S the upper (lower) row corresponds to 
normal (inverted) neutrino mass hierarchy. 

Besides 6*13 and S, from the global analysis of neutrino data we also recalculate the best fit values and ranges allowed 
for all the other neutrino oscillation parameters. Our results are summarized in Fig. [5] and Table [I] 

Comparing with our previous results we see that the inclusion of the new reactor and long-baseline data does not 
have a strong impact on the determination of the solar neutrino oscillation parameters, which are already pretty 
well determined by solar and KamLAND reactor data. The differences between the results in Table Q] and those 
in Table I in J3] are due to the different treatment of reactor data. Indeed, motivated by the so-called "reactor 
antineutrino anomaly" [ijj, old data from reactor experiments were included in the analysis in The dependence 
of the determination of solar neutrino oscillation parameters sin 2 6*12 and Am^ upon the details of the reactor data 
analysis has already been discussed in detail in Ref. [l|. 

Concerning atmospheric neutrino parameters, the best fit values for the atmospheric mass splitting parameter Am^ 
in Tab. U have been shifted to somewhat larger values compared to our previous results in Ref. Q. This is mainly due 
to the new MINOS disappearance data in Ref. [ljj, that prefer values for the mass splitting parameter larger than in 
their previous data release in |3o1 |. The precision in the determination of Am^ has also been improved thanks to the 
new long-baseline neutrino data. Thus, at 3c we find approximately a 8% accuracy in the determination of Ato^, 
while a 12% accuracy was obtained in Q at 3a. For the atmospheric mixing angle we note a slight rejection for 
maximal values of 823- In particular, our global fit shows a preference for the mixing angle in the second octant. This 
preference is very weak for the normal mass hierarchy case, where a local best fit point at sin 2 #23 = 0.427 appears 
with Ax 2 — 0.02, so that a symmetric A\ 2 profile can be seen at middle-top panel of Fig. [5] . For inverted mass 
ordering however, the profile is more asymmetric and a local minimum for 623 appears in the first octant only at 
A% 2 ~ 1.5. Maximal mixing, i.e. 823 = 7r/4 is disfavoured at ~ 90% C.L. for both hierarchies. While the preference 
for non-maximal values of the atmospheric mixing angle comes directly from the new MINOS data, the choice of a 
particular octant comes from the interplay of long-baseline, reactor and atmospheric neutrino data, as we will discuss 
in detail in the next section. 

IV. DISCUSSION 

We now discuss the impact of the new long-baseline data and the atmospheric neutrino analysis in the determination 
of the atmospheric mixing parameter 823- As already stated above the new disappearance data from MINOS show a 
preference for non- maximal values of 6*23. Due to the smallness of the associated matter effects in MINOS, these data 
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FIG. 3: Upper panels: contour regions with A\ 2 = 1, 4, 9 in the sin 2 #23 - sin 2 #13 plane from the analysis of long-baseline 
(MINOS and T2K) + solar + KamLAND data (left panel), long-baseline + solar + KamLAND + new Double Chooz, Daya 
Bay and RENO reactor data (middle panel) and the global combination (right panel) for normal hierarchy. Lower panels, same 
but for (inverted) neutrino mass hierarchy. 



are octant-symmetric and therefore say nothing about the octant of the atmospheric mixing angle #23 ■ However, the 
interplay with long-baseline neutrino appearance and reactor antineutrino data breaks the octant-degeneracy, leading 
to a small preference for values of O23 smaller than 7r/4. This is seen in the left panels of Fig. [3] where we have plotted 
the allowed regions in the sin 2 #23 - sin 2 #13 plane from the combination of long-baseline (MINOS and T2K) with solar 
+ KamLAND. These data samples prefer 6*23 values in the first octant and the same holds for the case when new 
reactor data are included, see middle panels in Fig. [3j Up to this point all statistical analysis are in agreement 4ll. |42 1 . 



Nevertheless, when we then include atmospheric data in the global analysis, differences in the determination of #23 
arise due to the differences in the analysis of atmospheric neutrino data. In fact, one sees that the effect of combining 
with the atmospheric neutrino data maintains the preference for non-maximal values of #23 but leads to a shift towards 
sin 2 823 in the second octant for both mass orderings, as seen in the right panels in Fig. [3] In contrast Refs. 4l|, |42| 
find #23 in the first octant for both spectra. Note however that the preference for a given octant in our analysis is 
still rather marginal, and 623 values in the first octant appear with A\ 2 = 0.02 and 1.5 for normal and inverse mass 
hierarchy respectively. 

As stated in Section IIII B[ for our global fit we use the official Super-Kamiokande analysis of atmospheric neutrino 
data in Ref. |3s1 |. performed within the one-mass scale approximation. This analysis shows a preference for maximal 
Q23 mixing, although a small deviation of 623 to the second octant is found in the case of inverse mass ordering. 
The most recent analysis of the Super-Kamiokande collaboration presented in Neutrino-2012 [43j and performed with 
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the inclusion of solar mass splitting corrections is in agreement with their previous results (obtained without these 
corrections) . In this case a small preference for # 2 3 m the first octant for the normal spectrum, and for second octant 
for the inverted one is found, with maximal mixing well inside the one sigma range. In contrast, the analyses of 
Refs. 



41 



42j, updated after Neutrino-2012, find a global preference for #23 in the first octant and exclude maximal 
mixing at the 2a level (for normal hierarchy), in qualitative agreement with each other, though the agreement is not 
perfect at the quantitative level. Both of these analyses are at odds with the latest Super-Kamiokande atmospheric 
neutrino data analysis in [4^ . At the moment it is not clear what is the origin of this discrepancy. 

The impact of the atmospheric neutrino analysis upon the determination of #23 is very visible, therefore in order 
to get a robust measurement of the atmospheric mixing angle it is crucial to clarify the origin of the discrepancies 
among the various analysis of atmospheric data. 



V. CONCLUSION AND OUTLOOK 

We have updated the global fit of neutrino oscillation parameters including the recent measurements of reactor 
antineutrino disappearance reported by the Double Chooz, Daya Bay and RENO experiments, as well as latest 
MINOS and T2K appearance and disappearance results, as presented at the Neutrino-2012 conference. We have 
found that the preferred global fit value of £> 13 is sin 2 13 = 0.0246(0.0250) for normal (inverted) neutrino mass 
hierarchy, while sin 2 #13 = is now excluded at 10.2a. There is reasonable agreement with the results of other global 
analyses [4l|, 0|, except for the atmospheric neutrino mixing parameter. We find that the global analysis pushes the 
atmospheric mixing angle sin 2 #23 best fit value towards the second octant for both neutrino mass orderings. This 
hint, however, is still quite marginal and first-octant values of #23 are well inside the la range for normal hierarchy 
and at 1.2cr for the inverted spectrum. Independent phenomenological analyses of atmospheric neutrino data in 
Refs. 41, 2] obtain a preference for mixing angle in the first octant for both mass hierarchies. Moreover, the new 
official Super-Kamiokande analysis in Ref. [43| with full three flavour effects gives a somewhat weaker preference for 
non-maximal #23 mixing, together with a correlation between the neutrino mass ordering and the preferred octant 
for 623- The origin of this discrepancy which crucially affects the determination of the atmospheric mixing angle is 
not yet clear. The impact of the new reactor and long-baseline accelerator measurements upon the solar neutrino 
oscillation parameters is completely marginal, the results are summarized in Table |H 

During the summer this year the Daya Bay Collaboration will complete the designed number of detectors by adding 
one detector in the far hall and other one in one of the near halls, re-starting the data taking after summer with eight 



neutrino detectors. After 3 years of operation the uncertainties on sin 2 2^13 will be reduced from 20% to 4-5% 44 1 
Needless to say that a good determination_of a sizeable 6*13 value will be a crucial ingredient towards a new era of CP 
violation searches in neutrino oscillations 
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46j and will also help determining the neutrino mass hierarchy. 



Note added: During the refereeing process of the present work, a new analysis has appeared in Ref. [47|. The authors 
obtain a preference for non-maximal values of #23 at the 1.7-2<r level as a result of the new MINOS disappearance 
data. Their new atmospheric neutrino data analysis results in a reduced sensitivity to the #23 octant in their global 



fit, in better agreement with the analysis by the Super-K collaboration than their previous results given in [42 1 



Acknowledgments 

This work was supported by the Spanish MINECO under grants FPA201 1-22975 and MULTIDARK CSD2009- 
00064 (Consolider-Ingenio 2010 Programme), by Prometeo/2009/091 (Generalitat Valenciana) and by the EU ITN 
UNILHC PITN-GA-2009-237920. M.T. acknowledges financial support from CSIC under the JAE-Doc programme, 



9 



co-funded by the European Social Fund. 



T. Schwetz, M. Tortola and J. W. F. Valle, New J. Phys. 13, 063004 (2011), [1103.0734]. 
T. Schwetz, M. Tortola and J. W. F. Valle, New J.Phys. 13, 109401 (2011), [1108.1376]. 
DOUBLE-CHOOZ Collaboration, Y. Abe et al, Phys.Rev.Lett. 108, 131801 (2012), [1112.6353]. 
DAYA-BAY Collaboration, F. An et al, Phys.Rev.Lett. 108, 171803 (2012), [1203.1669]. 
RENO collaboration, J. Ahn et al., 1204.0626. 

J. Schechter and J. W. F. Valle, Phys. Rev. D22, 2227 (1980); W. Rodejohann & J. Valle, Phys.Rev. D84, 073011 (2011) 
K. Nakamura et al, Journal of Physics G: Nuclear and Particle Physics 37, 075021 (2010). 



CHOOZ Collaboration, M. Apollonio et al., Eur. Phys. J. C27, 331 (2003), |hep-ex/0301017 



Palo Verde Collaboration, F. Boehm et al, Phys. Rev. D64, 112001 (2001), [hep-ex/0107009 



T2K Collaboration, K. Abe et al, Phys.Rev.Lett. 107, 041801 (2011), [1106.2822]. 

MINOS Collaboration, P. Adamson et al, Phys.Rev.Lett. 107, 181802 (2011), [1108.0015]. 

M. Gonzalez-Garcia, M. Maltoni and J. Salvado, JHEP 1004, 056 (2010), [1001.4524]. 

G. Fogli, E. Lisi, A. Marrone, A. Palazzo and A. Rotunno, Phys.Rev. D84, 053007 (2011), [1106.6028]. 

M. Ishitsuka's talk at Neutrino 2012 conference, Kyoto, June 2012 Y. Abe et al. larXiv: 1207.6632 1 [hep-ex] 

D. Dwyer's talk at Neutrino 2012 conference, Kyoto, June 2012 

R. Nichol's talk at Neutrino 2012 conference, Kyoto, June 2012 

K. Abe et al. [T2K Collaboration], Phys. Rev. D 85, 031103 (2012) arXiv:1201. 13861 [hep-ex]]. 
T. Nakaya's talk at Neutrino 2012 conference, Kyoto, June 2012. 



V. Kopeikin, L. Mikaelyan and V. Sinev, Phys.Atom.Nucl. 67, 1892 (2004), [hep-ph/0410100 
T. Mueller et al, Phys.Rev. C83, 054615 (2011), [1101.2663]. 
K. N. Abazajian et al, larXiv: 1204.53791 [hep-ph]. 



P. Vogel and J. F. Beacom,z Phys. Rev. D60, 053003 (1999), [hep-ph/9903554 . 
B. T. Cleveland et al, Astrophys. J. 496, 505 (1998). 

F. Kaether, W. Hampel, G. Heusser, J. Kiko and T. Kirsten, Phys.Lett. B685, 47 (2010), [1001.2731]. 
SAGE Collaboration, J. N. Abdurashitov et al, Phys. Rev. C80, 015807 (2009), [0901.2200]. 

G. Bellini et al. [ The Borexino Collaboration], Phys. Rev. Lett. 107 (2011) 141302 [arXiv: 1104. "18161 [hep-ex]]. 
Super-Kamkiokande Collaboration, J. Hosaka et al., Phys. Rev. D73, 112001 (2006), [hep-ex/0508053| . 

J. P. Cravens et al. [Super-Kamiokande Collaboration], Phys. Rev. D 78 (2008) 032002 [arXiv:0803 i43T2l ]hep-ex]]. 

Super-Kamiokande Collaboration, K. Abe et al, Phys.Rev. D83, 052010 (2011), [1010.0118]. 

SNO Collaboration, B. Aharmim et al, Phys. Rev. Lett. 101, 111301 (2008), [0806.0989]. 

SNO Collaboration, B. Aharmim et al, Phys. Rev. C81, 055504 (2010), [0910.2984]. 

A. Serenelli, S. Basu, J. W. Ferguson and M. Asplund, Astrophys. J. 705, L123 (2009), [0909.2668]. 

T. Schwetz, M. Tortola and J. W. F. Valle, New J. Phys. 10, 113011 (2008), [0808.2016]. 

The KamLAND Collaboration, A. Gando et al, Phys.Rev. D83, 052002 (2011), [1009.4771]. 

Super-Kamiokande Collaboration, R. Wendell et al., Phys.Rev. D81, 092004 (2010), [1002.3471]. 

MINOS Collaboration, P. Adamson et al., Phys.Rev.Lett. 106, 181801 (2011), [1103.0340]. 

MINOS collaboration, P. Adamson et al, Phys.Rev.Lett. 107, 021801 (2011), [1104.0344]. 

MINOS Collaboration, P. Adamson et al, Phys. Rev. Lett. 108, 191801 (2012), [1202.2772]. 

M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 6, 122 (2004), |hep-ph/0405172 

G. Mention et al, Phys. Rev. D83, 073006 (2011), [1101.2755]. 

G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo and A. M. Rotunno, Phys. Rev. D 86 (2012) 013012 
|arXiv: 1205.52541 [hep-ph]] . 

T. Schwetz's talk at "What's v" workshop, Florence, June 2012. 
Y. Itow's talk at Neutrino 2012 conference, Kyoto, June 2012. 

J. Cao's talk at nuTURN2012 - Neutrino at the Turning Point, Gran Sasso, May 2012. 

H. Nunokawa, S. J. Parke and J. W. F. Valle, Prog. Part. Nucl. Phys. 60, 338 (2008). 

ISS Physics Working Group, A. Bandyopadhyay et al, Rept.Prog.Phys. 72, 106201 (2009), [0710.4947]. 
M. C. Gonzalez-Garcia, M. Maltoni, J. Salvado and T. Schwetz. larXiv: 1209.30231 [hep-ph]. 



